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Abstract. Geometrical model of structure of the universe is 
examined to obtain analytical expression for the two points 
nonlinear correlation function. According to the model the ob- 
jects (galaxies) are concentrated into two types of structure 
1/-) elements - filaments and sheets. We considered the filaments 
ON ( similar to galaxy filaments ) simply as straight lines and the 
pH sheets ( similar to superclusters of galaxies ) simply as planes. 
£3 The homogeneously distributed objects are also taken into con- 
5 sideration. The spatial distribution of lines, planes and points is 
"^j" uncorrected. 

The nonlinear correlation function depends on four parame- 
ters and is similar to the observed and simulated ones for differ- 
CN ent samples. It describes quite well the correlation of galaxies, 
^ clusters of galaxies and dark matter distribution. 
q Possible interpretation of the parameters of nonlinear cor- 
l/-) relation function is discussed. 

On 

Key words: cosmology: theory - large scale structure of the 
O ^niverse - correlation function 
I 

o 

& 

1. Introduction 

The two point correlation function £(| n - r 2 |) = £(r) is 
the most popular tool for the investigation of spatial object dis- 
tribution both in the observational and simulated catalogues. 
The great potential of this method is known and it can be illus- 
trated by the new interesting results obtained at the last years. 
Thus, Jorgensen et al. (1993, 1994) shown that oscillations in 
the spectra of primordial perturbations give rise to some inter- 
esting peculiarities in the correlation function. Let us remind 
also the paper of Mo et al., (1992) where characteristic scales 
were found in observational catalogues. Recently Peacock & 
Dodds (1994) reconstructed the initial power spectrum using 
the correlation analysis of different samples of galaxies and 
galaxy clusters. This list can be continued. 
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The main drawback of correlation method is, of course, its 
phenomenological character. Indeed, the interpretation of non- 
linear correlation function, £ > 1, is ambiguous. The same 
correlation function can be generated by different spatial point 
distributions. For example, the similar correlation functions are 
generated by Soneira-Peebles (1978) model, by Voronoi tessel- 
lation (van de Weygaert, 1991), by the 'shell' model (Bahcall 
et al., 1989) and by many simulations with HDM, CDM and 
other power spectra. It is important that the spatial matter dis- 
tribution in these models is very different. However the shape 
of nonlinear correlation functions is the same for a wide class 
of power spectra (Demianski & Doroshkevich (1992)). 

The main properties of observed correlation function still 
have no physical explanation. Thus, for all observational sam- 
ples the shape of the correlation function (for galaxies and 
clusters of galaxies) is approximately the same: 

t{r) = (r/r c )"^ (1.1) 

with the power index 7 ss 1.8. However, the correlation radius, 
r c , as a rule, depends on the sample. Thus, for IRAS, CfA and 
APM catalogues r c ss (5 - 6) h~ x Mpc (Davis et al., 1988, 
Loveday et al., 1992). But for the sample of dwarf galaxies 
(Thuan et al., 1991) r c « 4 h~ l Mpc and r c « Ihh^Mpc for 
the Great Wall region (Ramella, Heller & Huchra, 1992). Of 
course, these variations are produced by strong peculiarities of 
the samples. However, adequate physical explanation of these 
phenomena are still lacking. 

Let us also note that the strong correlation of clusters of 
galaxies has no physical explanation. Of course, similar corre- 
lation has been simulated (see, e.g., Klypin & Rhee, 1994, Croft 
& Efstathiou, 1994) and, therefore, the parameters of correla- 
tion function can be caused by the power spectrum. However, 
at present there is a strong limit to what can be done in this area. 

There are different interpretations of the correlation func- 
tion. Thus, it can be considered as the evidence in favor of 
hierarchical clustering (Soneira & Peebles, 1978) or in favor of 
fractal character of galaxy distribution (Coleman & Pietronero, 
1992). Here we interpret the shape of correlation function as the 
evidence in favor of strong concentration of galaxies in struc- 
ture elements. Evidently, additional information is necessary 
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for a proper interpretation of the correlation function. Never- 
theless, the correlation analysis gives the important quantitative 
information about the spatial matter distribution and it cannot 
be ignored in theoretical considerations. 

The simplest model of the nonlinear correlation function 
is provided by a geometrical model. According to this model 
the nonlinear correlation function is generated by the concen- 
tration of significant fraction of the objects along lines and 
planes. The great potential of this approach was emphasized 
by Zeldovich (see, e.g. Shandarin et al. 1983). This model is 
suitable to describe the Large and Super Large Scale Structure 
(LSS and SLSS) in the observed distribution of galaxies in deep 
galaxy surveys and the dark matter (DM) distribution in dynam- 
ical N-body simulations (Buryak, Doroshkevich & Fong, 1994, 
hereafter BDF, Doroshkevich et al., 1995, hereafter DFGMM, 
Doroshkevich, 1995). In the following we use this model to 
construct the nonlinear correlation function. 

In the geometrical model filaments are considered simply 
as lines and sheets simply as planes. To be more specific, we as- 
sume that a filament is identified with a line segment and a sheet 
corresponds to a finite planar region. If filaments and sheets are 
large enough we neglect the edge effects in our analysis. 

This model can be employed to describe the spatial dis- 
tribution of various objects. Thus, for the observed catalogues 
a 'filament' is represented by a chain of galaxies or clusters 
of galaxy, a 'sheet' is represented by superclusters or richest 
filaments of galaxies. For the DM distribution in simulated cat- 
alogues a 'filament' can be identified with strongly extended 
over dense region, whereas a 'sheet' with high dense wall at the 
boundary of voids or under dense regions. 

This approach is not general and it has its limitations. It 
is an intermediate step between the local description of matter 
distribution by means of density and velocity fields and such 
global methods as the percolation or fractal analyses. The ge- 
ometrical model works at scales comparable to the correlation 
radius and, therefore, it can be used to describe the nonlinear 
correlation function. 

Random distribution of straight lines can be characterized 
by the surface density, er/, i.e., the mean number of lines inter- 
secting a unit area of arbitrary orientation. Similarly, random 
distribution of sheets can be characterized by the linear density, 
a s , i.e., the mean number density of sheets crossing an arbitrary 
straight line. Of course, we have to include also some fraction 
of homogeneously distributed points. It is characterized by the 
3D number density, n^. Parameters of the correlation function 
- the power index, 7, the correlation radius, r c , and the zero 
point, ro, - are related to structure parameters. 

In the geometrical model the three and four point correlation 
functions can be analytically calculated as well. However, these 
expressions are not so simple and, actually, the main advantage 
of analytical approach is lost. Simple numerical simulations of 
such matter distribution seem to be more perspective ways to 
examine the three and four point correlation functions. 

This paper is organized as follows. In section 2 the model 
of the correlation function is formulated. In Section 3 we exam- 
ine this model and obtain the two point correlation function of 



galaxy distribution. In Section 4 we reconstruct the DM correla- 
tion function using some information about the DM distribution 
from dynamical N-body simulations. In Section 5 we propose 
a set of models for the cluster-cluster correlation function. In 
Section 6 and 7 the evolution of correlation function and our 
results are discussed. 

2. The analytical model of correlation function 

In this section we consider the correlation function for three 
simple objects distributions: a homogeneous point distribution 
and randomly distributed straight lines and planes. Further we 
construct the composite correlation function including all three 
types of point distribution. 

2.1. Correlation function of homogeneous point distribution 

The simplest model of the object distribution is a homogeneous 
Poissonian distribution with a mean number density < n g >= 
n h . As it is known, in this case the number of actual pairs, 
DD, at the distance (r, r + dr) is identical to the number of 
Poissonian pairs, PP, at the same distance 

DD = PP = <n g > 4irr 2 dr = Aim h r 2 dr (2.1) 

Therefore, in this model we obtain the trivial result: 

ew = (2.2) 

2.2. Correlation function of randomly distributed straight lines 

The more interesting model is a system of randomly distributed 
identical straight lines. It is characterized by the surface density 
of lines, er/. We assume that objects are distributed homoge- 
neously along lines with a linear number density rii . 

Let us consider a sphere with a radius r which has an ob- 
ject and, therefore, a straight line at its center . In thus case 
the main characteristics of the point distribution are: the mean 
number of lines, Nf, intersecting the sphere, the mean length 
of intersection, If, and the mean 3D number density of objects, 
< n g > inside the sphere. These parameters can be found by 
the standard method (see, e.g., BDF). Here we present only the 
final expressions: 

Nf = 0.5(7/5 = 2 7 rr 2 <7 / , l f = AV/S = \r 

<n>= mlfNf/V = 2(7/711 (2.3) 

where V and S are the volume and surface area of the sphere. 

However, it must be taken into account that a line is crossing 
through the center of sphere. Thus, we have to use the condi- 
tional mean number of lines intersecting the sphere. This value, 
iV/, can be written in the following form: 

1 00 N 

«; = ?T 3!wg«'<+.W = a3k)- 1 <2 - 4) 
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where Wk(Nf) is a probability to find k lines intersect- 
ing the sphere. Therefore, using the Poissonian probability 



N 



L e N > we obtain: 

4 

DD = 2n L dr + -mdrN), 

O 

PP = Stirr 2 o f7i\dr = 4Nfri\dr, 



Therefore, for the composite model we obtain 



DD = fh ( Anr 2 drnh + —n\drNj + -irrdrn^ N ( 



6N f 



e N f - 1 



(2.5) 



Thus, for small r (TV/ « 1) f = (27V/)- 1 oc r~ 2 , and for 
large r (N f » 1) f (r) ss -2/3. 

2.3. Correlation function of randomly distributed sheets. 

This case can be treated in the same way as a previous one. 
We assume that the objects are distributed homogeneously on 
identical sheets with a surface density n 2 . 

Let us consider again a sphere with radius r which has in 
its center an object and, therefore, a sheet. For this model the 
mean number of sheets, N s , intersecting the sphere, the mean 
area of intersection, S s , and the mean 3D number density of 
objects inside the sphere, < n g >, are: 

2 

N s = 4cr s r, S s = -Trr , < n g > = 2a s n 2 (2.6) 

The conditional mean number of sheets intersecting the sphere, 

TV* is 



+ // (^2n\dr + + -riidrNJ + -Trrdrn 2 N s + 4Trr 2 drnf^j 

( 4 4 

+ f s I 2irrdrn2 + —n\drNf + —irrdrn^Nl + 4irr 2 drnh 

PP = 4irr 2 dr < n >= 4irr 2 drn h / f h 
= 4Trr 2 drrif / // = 4Trr 2 drn s / f s 

where it is assumed that a) the sphere is centered on an ob- 
ject from the Poisson population, and b) the sphere is centered 
on an object from filament and sheet populations. Finally, the 
analytical expression for the correlation function is 



6N f 



e N f -1 



' 1 + 37V ' 



-(/ s + 2//)/3 

A + // + /. = i 



(2.10) 



Let us remind that according to (2.3) and (2.6) Nf = 2irr 2 af, 
N s = 4a s r and, therefore, for smallr £(r) oc r -2 . This com- 
posite model reduces to (2.2), (2.5) and (2.8) when = 1, 
// = 1 and f s = 1 correspondingly. 



1 °° TV 
^ s = (1 — W ) ^ (-^*) = (i - w ) ~ ^ Some properties of the correlation function. 



Therefore, in this model 



DD = 2-Krdrri2 + —TrrdrrizN*, 
PP = %'Kr 2 dra s n2 = 4N s ri2Trrdr, 



2N„ 



1 



(2.8) 



and, therefore, ^ = 1/N„ oc r 1 for small r (N s « 1), and 
£(r) « -1/3 for large r (N s » 1) 

2.4. The composite model of correlation function. 

It is clear that the simple models considered above cannot be 
used to describe any realistic catalogues. It is more reasonable 
to use composite model in which all three types of structure 
elements (sheet, filament and homogeneous populations) are 
mixed with the fraction factors f s , // and correspondingly. 

The mean density of particles and fraction factors can be 
written as 

< n >= n h + 2a jrii + 2a s n 2 = n h + nj + n s (2.9) 
fh=n h /<n>, ff = 2a f nx/ <n >, f s = 2a s n 2 /<n> 



Let us consider in more detail the relationship (2.10) and 
find the main properties of correlation function, £(r). Natu- 
rally, our model cannot be applied for large distances, r » 
r c , £(r) << 1, where the matter distribution 'remembers' the 
initial perturbations and, therefore, the small correlations in the 
distribution of structure elements are essential. Let us remind 
that we assumed the uncorrelated spatial distribution of struc- 
ture elements and, so, the correlation function (2.10) describes 
only the matter concentration into structure elements. 

Thus, we will concentrate our attention on the nonlinear 
region where the correlation function is large enough and our 
approach is justified. As a rule, in this region, for the analytical 
consideration we can use the simpler version of (2.10): 



£(r) = (r//r) 2 + 2(r./r) - (2.11) 

r 2 = f 2 f/4ira f , r s = f 2 /8a s , ^ = (/. + 2//)/3 

where we assumed that Nf < 1, N s < 1. 

The correlation radius, r c , is defined by the relation 



1, 



(1 + Coo) 



(i + V^Te) , (2.i2) 



C 2 = (l + £ooXr//r s ) 2 
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Another important parameter of the correlation function, 
the zero point, ro, can be found from 

f(ro) = 0, „ = i(, + ^ 1 + _^_ C .) M 

For the important case r s > 77 , £ <C 1 these expressions 
simplify to 

r c « 2r s /(l + £oo), r « 2r s /£ 00 

r /r c - 1 « l/^oo (2.14) 

It establishes the interesting relations between the main pa- 
rameters of the correlation function and characteristics of the 
structure. 

For the model (2.11) near the point r = r c the power index 
of the correlation function 7 is a weak function of r: 

7 = 7c(l + &(--!) + ..•) (2-16) 



7c = 2(l + Coo) ^—- (2 - , l + ^oo<7c<2(l + Coo) 

& = V7c - 1 - £00 (3 - 7c " 2(1 + £oo)/7c) 

It is interesting that in two limiting cases, ( <C 1 and ( > 1 
the power index is connected with r c and r as follows: 

7c «(l - r c /r )-\ C«l (2.17a) 

7c « 2(1 - r 2 Jr*)-\ C»l (2.176) 

3. The correlation function of galaxy distribution 

In this section we consider the composite model (2.10) with 
the structure parameters found from the observed galaxy distri- 
bution. We use the estimates of the main structure parameters 
obtained by BDF from deep galaxy surveys. 

Thus, according to BDF, the distribution of galaxy filaments 
can be characterized by the surface density 

07 « 0.5 x W- 2 h 2 Mpc- 2 (3.1) 

In this case the SLSS elements can be considered as the sheet 
population. The characteristic scale of the SLSS observed from 
two different types of the surveys gives: 

a s « 2 x lO- 2 hM P c-\ a s « 0.77 x lO^hMpc' 1 . (3.2) 

To test the possible shape of the composite correlation func- 
tion for these structure parameters two models with the same 
value of o f and different parameters a s , ff and f s were cal- 
culated. For both models we assumed that all galaxies are con- 
centrated into the structure elements and, thus, f^ = 0. 

The functions £(r) are shown in the Fig. 1 forthe correlation 
radius r c = 5.5h~ 1 Mpc. The parameters of the models are 



^^^^ 1 1 1 1 1 1 








1 1 



1. 10 

r fa -1 Mpc 



Fig. 1. Correlation function (2.10) for galaxies. Parameters are pre- 
sented in Table I. Straight lines correspond to the function (3.3). 

presented in Table I. The straight line in Fig. 1 corresponds to 
the function 

«r) = (rc/r) 1 - 8 (3.3) 

These results illustrate the potential of the composite model and 
allow us to evaluate the galaxy fractions f s and ff in different 
elements of the structure. 

Table I 



CT S 10- 2 


/. 


n 2 h 


cr/10- 2 


// 


ri\h 


roh- 1 


hMpc' 1 




Mpc' 2 


h 2 Mpc~ 2 




Mpc' 1 


Mpc 


2. 


0.84 


0.21 


0.50 


0.16 


0.16 


15.2 


1. 


0.58 


0.29 


0.50 


0.42 


0.42 


14.8 


1. 


0.56 


0.28 


0.25 


0.44 


0.88 


13.8 



There is however a disagreement between estimates (3.1), 
(3.2) and assumptions of the model (2.10). The value of the 
filament surface density, Of, used above was obtained by BDF 
as a low limit of size of the structure when the poorest fila- 
ments were taken into account. However, in the model (2.10) 
all filaments are considered as identical and, therefore, the sur- 
face density of the 'average' filaments should be used. It could 
be smaller by a factor of 1.5 - 2 than the accepted value. To 
test such 'correction' in Fig. 1 the correlation function for the 
smaller surface density Of ( the model 3 in Table I ) is shown as 
well. In this case the correlation radius r c = t z>.^>hr 1 Mpc but 
the fractions f s and ff have different values than above. 

Other values of the correlation radius r c can be obtained for 
other values of the fractions ff, f s . We conclude that, actually, 
the correlation radius depends on many structure parameters in 
the investigated catalogue. 
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The estimates of the mean linear density of galaxies in 
filaments, n\, and surface density of galaxies in sheets, are 
also presented in Table I for < n g >« lO~ 2 h 3 Mpc~ 3 . The 
mean separation of galaxies in the structure elements, 1 and 
n^ 1 ^ 2 are less than < n > -1 / 3 by a factor of 2-5. 

4. The correlation function of dark matter distribution. 

For simulated DM distribution, as a rule, the correlation func- 
tion of DM spatial distribution is characterized by the same 
power index, 7 « 1.8, as for galaxies and the correlation ra- 
dius, r c « 3.h~ 1 Mpc. Usually this difference in correlation 
radius is attributed to the 'biasing' factor, b. Analysis of LSS 
and SLSS in numerical simulations ( DFGMM) shows however 
that the spatial distribution of 'galaxies' and DM is very differ- 
ent. Therefore, it is interesting to test consistency between the 
parameters of the correlation function and DM structure on the 
basis of model (2.10). 

As it was found by DFGMM the DM distribution is dom- 
inated by the sheet like structure. Filamentary structure exists 
mainly in the smallest scale and it includes only a small fraction 
of the DM population. But, unlike in the galaxy distribution, the 
fraction of homogeneously distributed DM particles could be 
significant (10-20%). 

Thus, properties and parameters of the LSS for 'galaxy' 
and DM populations are very different. However, the shape 
and the power index of the correlation function are the same 
for both populations. The model (2.10) allows to obtain the 
correlation function with the reasonable correlation radius, r c ss 
3.h~ 1 Mpc, and with the power index, 7 « 1.8 for the range 
parameters presented in Table II. 

Table II 



t s io- 2 


/. 


ff/ io- 2 


// 


Teh' 1 


roh- 1 


hMpc' 1 




h 2 Mpc~ 2 




Mpc 


Mpc 


4. 


0.88 


1.0 


0.12 


3.0 


8.5 


4. 


0.39 


0.25 


0.61 


3.0 


5.2 


3. 


0.72 


0.25 


0.03 


3.0 


11.0 



Estimates of a s and Of found by DFGMM were used. The 
corresponding correlation functions are shown in Fig. 2. 
The straight line in Fig. 2 corresponds to the function 

WO = (r c /r) 18 (4.1) 

withr c = Z.hT l Mpc. 

Notice again that in the model (2.10) all filaments and 
sheets were assumed to be identical while the direct analysis 
(DFGMM) shows distribution of the structure elements with 
richness. Thus, the model (2.10) should be also corrected for 
this factor. However, in this paper we are interested in depen- 
dence of the correlation function on other structure parameters. 
For example, the model 2 in Fig. 2 and in Table II corresponds 
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Fig. 2. Correlation function (2. 10) for DM distribution. Parameters are 
presented in Table II. Straight lines correspond to the function (4.1). 

to the smaller surface density of DM filaments. The model 3 
in Fig. 2 and in Table II illustrates the influence of the ho- 
mogeneous background of DM particles (f h = 0.25) on the 
correlation function. 

This sample of correlation functions demonstrates that the 
same correlation function can be generated by different DM 
distribution. 

5. The correlation function of clusters of galaxies. 

The observed cluster-cluster correlation function is character- 
ized by the power index, 7, by the correlation radius, r c , and by 
the zero point, ro : 

7 « 1.8, r c « 2Qh~ 1 Mpc, r « (50 - 60)/i _1 Mpc 

(5.1) 

(see, e.g.,Daltonetal., 1992, Peacock & West, 1992, Scaramella 
et al., 1993). The mean number density of clusters, n c i, and the 
mean separation of the clusters, D c i, are 

n cl « 2.4 x 10" 5 /i 3 Mpc -3 , 

D cl = n d 1,Z « 35/T 1 Mpc. (5.2) 

The values of r c , ro, n c ; and D c i strongly depend on the 
sample. There are also other estimates (Bahcall, 1988, Postman, 
Huchra & Geller, 1992) 

n cl « 6 x 10" 6 /i 3 Mpc -3 , D cl « 55/i _1 Mpc. (5.3) 

Estimates of the correlation radius, r c , and the zero point, ro, for 
some samples may be as high as 50h~ 1 Mpc and 100h~ 1 Mpc 
correspondingly (Bahcall, 1988). 

Table III 



6 



O. Buryak, A. Doroshkevich,: Correlation function as a measure of the structure 



rr If!" 2 

JiMpc^ 1 


Is 


n„i n — 2 

h 2 Mpc^ 2 


rr <-1 n — 4 

/ 

h 2 Mpc^ 2 


f , 
Jf 


Ti\ h 
Mpc^ 1 


Mpc 


n ?n 


n 3d 


n 94 

U.Z4 


n R4 

U.04 


n zti 


U.Uo 


f;i 
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n ^4 

U.04 
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4fi 
40. 


0.20 


0.36 


0.21 


0.16 


0.20 


0.15 


64. 


0.15 


0.27 


0.22 


0.09 


0.17 


0.23 


68. 


0.10 


0.27 


0.32 


0.50 


0.33 


0.08 


60. 



Therefore, it can be expected that the spatial distribution of 
clusters can be described by our analytical model. It is signifi- 
cant that a numerical 'shell' model (Bahcall et al, 1989) which is 
close to the model (2. 10) in many respects is consistent with the 
observed cluster correlation. Moreover, the observed parame- 
ters r c , r and 7 (5.1) approximately satisfy the relationship 
(2.17a). 

We present here five models of correlation function. Their 
parameters are listed in Table III. Fig. 3 shows the corresponding 
correlation functions. For all models the correlation radius is 
equal to 2Qh~ 1 Mpc, and we used data (5.2) to obtain the mean 
number density of clusters in the structure elements. These 
models show that, indeed, it is possible to interpret the observed 
and simulated cluster-cluster correlation functions as a result of 
the cluster concentrations into structure elements of the SLSS. 

Of course, the spatial distribution of clusters should be in- 
vestigated in more detail to reveal the quantitative characteris- 
tics of their distribution. However, the small values of of and 
a s in Table III show that it is necessary to use very large cata- 
logues in order to get statistically significant estimates. On the 
other hand, found here mean separation of clusters in sheets 
D s cl « (20 - 25)h~ 1 Mpc is close to that in the homogeneous 
population 

D h cl « 35/~ 1/ V 1 Mpc « Ahh^Mpc. (5.4) 

This means that it is very difficult to identify sheets in obser- 
vational catalogues. For the filament population the expected 
mean separation of clusters is close to hhr l Mpc and it is more 
easier to identify it. 

6. Evolution of the correlation function. 

Evolution of the correlation function has been simulated for 
different spectra of the primordial perturbations and, in general, 
it can be described as a growth of amplitude (or correlation 
radius, r c ) and power index, 7. On the other hand, the direct 
analysis of the DM structure evolution (DFGMM) shows that 
there are three main evolutionary periods: 

1 . The structure arises during a short period which is defined 
by the amplitude of density perturbations. At the same time the 
nonlinear part of correlation function is formed. 

2. Later the typical scale of the LSS grows slowly. 

3. During the last period of evolution sheets are transformed 
to filaments and further to a system of knots. 

Therefore, as a rule, // is an increasing functions of time 
while f s is decreasing one. Expressions (2.12), (2.16) show 
that this picture of structure evolution is in a qualitative agree- 
ment with evolution of the correlation function. However, the 
presently available data does not allow us to test this evolution- 
ary history more quantitatively. 

For the cluster - cluster correlation function the evolution 
is very slow because the mean velocity of clusters is relatively 
small. Probably, the continuous formation of new clusters can 
be considered as an essential factor in the evolution of the spatial 
distribution of clusters. 




10 100 
r IT 1 Mpc 

Fig. 3. Correlation function for the clusters of galaxies. Parameters 
are presented in Table III. Straight lines correspond to the function 
£ci = (r c /r) 18 withr c = 20/T 1 Mpc. 

Significant fraction of clusters forms chains in superclusters 
of galaxies which are similar to the well known chain in the 
Perseus supercluster (rich clusters Abel 426, 347 & 262). An 
example of clusters walls has been demonstrated by Kopylov 
et al., (1988). A catalog of clusters of clusters of galaxies was 
prepared by Bahcall & Soneira (1984). The linear size of such 
clusters of clusters could be as large as lOO/i -1 Mpc. 

On the other hand, the percolation analysis gives the value 
P « 3 (Klypin, 1988) for the percolation parameter which is 
close to the Poissonian one. This fact shows that the cluster pop- 
ulation has a less pronounced structure and there is a significant 
fraction of isolated clusters. 
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7. Discussion. 

The main aim of this paper is to show that the nonlinear part 
of observed and simulated correlation functions can be gener- 
ated by the concentration of the objects (galaxies, DM particles 
and clusters of galaxies) into structure elements - filaments and 
sheets. Of course, the converse is not true and various spa- 
tial distributions of objects can generate the same correlation 
function. It is essential that in order to construct the model of ob- 
served correlation function we have to use no less than two types 
of structure elements (filaments and sheets) and, therefore, no 
less than three structure parameters. Probably, the nonhomo- 
geneous matter distribution into structure elements produces 
similar effects in numerical models (Bahcall et al., 1989, van 
de Weygaert, 1991). 

Thus, we can conclude that the parameters of the nonlin- 
ear correlation function depend on many structure parameters 
and, therefore, the interpretation of the correlation function is 
ambiguous. 

The model (2.10) gives a good approximation of the ob- 
served correlation function in nonlinear area for a wide range 
of structure parameters. Moreover, using the observed values of 
basic structure parameters, such as the surface density of fila- 
ments, af, and linear density of sheets, a s , we can also evaluate 
the fraction of objects concentrated in different types of the 
structure elements and find the reasonable mean separation of 
objects in structure elements. These results show that in many 
instances the model (2.10) corresponds to the actual situation. 

The analysis of Sec. 3 & 4 shows that the biasing parame- 
ter, b, describing the difference between the spatial distribution 
of DM and 'galaxies' actually depends on many structure pa- 
rameters. Therefore, the biasing parameter is not an universal 
characteristic of matter distribution (or DM model), but, rather, 
it is a characteristic of the formation of galaxies and structure. 
Thus, for galaxies and clusters the structure parameters have 
different origin and there is no reasons to expect the same bi- 
asing parameter for these objects. This conclusion is consistent 
with results of Peacock & Dodds (1994). Moreover, there are 
good reasons to think that the biasing parameter depends on the 
considered sample of objects and, so, it is different for over and 
under dense regions. 

The model (2.10) is very simple and there is a limit to what 
can be done in its framework. Thus, this model cannot explain 
the remarkably stable value of the power index 7 = 1.8 in 
all observed samples for all classes of objects. Of course, the 
DM and galaxy distributions evolve and therefore the value of 
7 characterizes the evolutionary period. However, evolution of 
the cluster distribution is slow (Croft & Efstathiou, 1994) and, 
perhaps, it manifests in the continuous cluster formation. Thus, 
in this case the power index characterizes rather the spatial 
distribution of regions where clusters form and emergence of 
the same value of 7 is surprising. 

As noted in the introduction, for some observed samples of 
galaxies the correlation radii are different although the power in- 
dex is the same. On the basis of the model (2.10) this fact can be 
accounted for by changing the structure parameters while main- 



taining the same (or close) fractions f s and ff. Probably, the 
universal dimensionless correlation function (Bahcall, 1988) 
can be interpreted by the same way. However, such explanation 
is not fully satisfactory and warrants further investigation. 

Sometimes the power law dependence of the correlation 
function (1.1) and the surprisingly universal value of 7 « 1.8 
are considered as the evidence in favor of fractal and multi- 
fractal properties of the matter distribution (see, e.g., Coleman 
& Pietronero, 1992). Actually, such approach may be effec- 
tively used to obtain additional quantitative characteristics of 
the object distribution. However, the previous discussion shows 
that these facts are not a good basis for the modification of 
fundamental concepts of cosmology. Various evolutionary his- 
tories can result in similar correlation and fractal properties and 
therefore interpretation of these properties is ambiguous. 

Traditionally, the strong correlation of clusters is described 
by the theory of linear Gaussian fluctuations in the smoothed 
density field (see, e.g., Holtzman & Primak, 1992, Croft & 
Efstathiou, 1994). However, such approach does not provide 
physical explanation of the strong cluster-cluster correlation 
and, so, it has rather phenomenological character. 

However, the results of Sec. 5 show that the observed cor- 
relation function can be generated by inclusion of the clusters 
into structure elements mean separation of which exceeds the 
correlation radius by a factor of 5 - 10. In this way we have to 
consider range of scales which are typical for the spatial dis- 
tribution of gravitational potential of perturbations and, thus, 
the strong correlation can be attributed to the concentration of 
clusters into elements of the SLSS. 

Concentration of clusters in the structure elements is con- 
sistent with general ideas of the nonlinear theory of gravita- 
tional instability. Indeed, as it has been argued by Buryak et 
al. (1992) and Kofman et al. (1992), the spatial distribution of 
SLSS elements is determined by the spatial distribution of the 
gravitational (or velocity) potential. It is therefore reasonable to 
expect that the spatial distribution of rich clusters is also deter- 
mined by the gravitational potential of primordial perturbations. 
Rich clusters should be concentrated mainly into a system of 
2D sheet and ID filament inside deep 'wells' or, partly, around 
very deep minima of potential. If this assumption is correct then 
the SLSS and the structure in the spatial distribution of clusters 
are generated by the same causes. In this case the strong clus- 
ter - cluster correlation is simply a consequence of preferential 
cluster formation in the SLSS elements. To test this hypothe- 
sis we have to study the cross correlation of the potential and 
clusters distribution in numerical simulations. Such an analy- 
sis can also determine the fraction of clusters distributed quasi 
homogeneously. 

However, although the concentration of clusters and galax- 
ies into structure elements is similar, what was emphasized, for 
example, by the universal dimensionless correlation function 
(Bahcall, 1988 ), this analogy has rather formal character be- 
cause the physical causes of the LSS and SLSS formation and 
their evolution are different. 
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